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ABSTRACT 


An alternating gradient magnetometer has been interfaced to 
a computer for the automation of data taking. Using a fast 
Fourier transform analysis system data can be acquired and 
processed in real time. Data is stored on disk and can be 
recalled for plotting and further analysis. With the 
addition of a simple liquid nitrogen cryostat, magnetization 
measurements can be carried out in the range from 300 K to 
77K. Results are reported on three different types of 
piezoelectric transducers. 
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Introduction 


The discovery of the ceramic superconductor Y-Ba-Cu-0 1 and 
its related compounds has led to a renewed interest in the 
general field of superconductivity. The ceramic 

superconductors, which have transition temperatures ranging 
from 90K to 120K and perhaps higher 2-4 have given rise to 
expectations of great technological utility. Unfortunately, 
most of these materials have low critical currents, which 
limits their current carrying abilities. This is a serious 
drawback for applications involving superconducting magnets 
and power transmission and storage. 

In order to understand the limitations to the critical 
current it is important to be able to characterize the 
magne- ir^ properties of a superconductor, as the two are 
intimately related . Measurements of magnetization and 
magnetic hysteresis allow one to determine the extent of 
.netic flux pinning, which is related to the microscopic 
critical current. 

'.n the current research an alternating gradient 
magnetometer, which was constructed in the Summer of 1988, 
was interfaced with a computer in order to automate the 
taking of magnetization data. 
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Objectives 


The objectives of this experiment were to interface the 
alternating gradient magnetometer to a computer and to write 
the program which would automate data taking. The program 
must be capable of reading the magnetization, temperature, 
and applied magnetic field, in real time, and storing them 
on a disk for later analysis. 
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Theory 


The internal magnetic field in a non-ferromagnetic material 
is described by the equation 


B = H + 4rrM 


( cgs units ) 


where B is the magnetic induction, H is the external applied 
new and M L the magnetization For homogeneous 
paramagnetic and diamagnetic materials M is _ a linear 

Eu^c trior of H. Thus, in a homogeneous material M - XH, whe.- 
X is the magnetic susceptibility. In a paramagnetio material 
X is positive and in a diamagnetic material it is negati_ . 
Superconductors are diamagnetic when in the superconducting 
state and (usually) paramagnetic or non-magnetic ( X ) - 

the normal state. So a graph of M as a .unction of Hsh 

t , dfraicht line of slope X. In a superconductor, this 
\ - i 1 1 be true as long as the applied field does not approacn 
the low r critical field Hc^ Once H C]L is exceeded, magnetic 
flux bee : ns to penetrate the sample, in the foCTti 
quantized flux lines. These flux lines become pinned m the 
camSIe, reducing the diamagnetic magnetization. The presence 
o- "pinned flux therefore, produces a hysteresis m the 
mignetTzation vs field curve.%he extent of the hysteresrs 
is 9 an indication of the extent of flux pinning. If an 
electric current is introduced in the sample, once i 
exceeds some critical value the flux lines will begin to 
move, giving rise to an electrical resistance. So, the 
extent of the hysteresis is also related to the critical 
currant, at which electrical resistance begins to appear. 

If a paramagnetic or diamagnetic material, is placed in an 
external, non-uniform, magnetic field, it experiences a 
force in the direction of the field gradient dH/dx. This 
force is directly proportional to the field gradient, and 
M, the magnetization. If one can measure this force then, 
knowing the field gradient, the magnetization can be 
computed . 

The theory and construction of an alternating gradient 
magnetometer has been described by Flanders and by the 
author ' . In the original design, the gradient field was 
swept sinusoidally and the force on the sample was measured 
usinq a piezoelectric transducer connected to a lock in 
amplifier. In its present form the piezoelectric output is 
detected by a fast Fourier transform analyzer. 
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Apparatus 


Magnetometer 


A block diagram of the magnetometer and associated circuitry 
is shown in figure 1. The major components will be discussed 
separately in the following paragraphs. A complete 
description of the field coils and vane can be found in 
reference 7 . 


Wave Analyzer 

The output of the piezoelectric transducer is processes by 
an Ono Sokki Mini FFT Analysis System. The FFT system is a 
small dedicated computer which performs a fast Fourier 
transform on the input signal and displays the wave spectrum 
n ® ar time. The amplitude of the peak corresponding 

to the driving frequency can be selected thereby eliminating 
other peaks caused by noise and overtones. The amplitude of 
the selected peak is read by the computer and stored in an 
array . 


Function Generator 

The driving current for the sweep coils is obtained from a 
Wavetek model 270 Function Generator. The drive voltage and 
frequency are input at the start of the program and written 
to the Wavetek by the computer. The frequency is set at the 
resonant frequency of the reed, or one of its overtones. The 
drive voltage is set to produce an rms current of 100 mA 
through the sweep coils . 


Current Source 

The current for the D.C. magnet is obtained from a Hewlett- 
Packard model 228A Current Source. The magnet current is 
ramped m steps of 0.10 Amp preceding each measurement of 
the signal amplitude. It was necessary to add a time delay 
of 2 to 4 seconds in order to allow the system to settle 
down between measurements. A further provision was added to 
allow the sample to be cooled initially in a preset magnetic 
field, in order to analyze the extent of flux pinning. 
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Figure 1. Magnetometer Block Diagram. 
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Computer 

The HP 9133 computer was programmed in Basic to control the 
magnetometer and collect data on the magnetization, magnetic 
field, and temperature, through the GPIB bus. Data could be 
saved to disk and printed out or plotted. An HP 631 printer 
and an HP 7470A plotter were connected to the computer lor 
these purposes . 


Temperature Measurements 

The temperature sensor is a Copper-Constantan thermocouple 
which is glued to the end of the vane, next to the sample, 
with GE 7031 varnish. The thermocouple output goes to an HP 
Model 3451A programmable multimeter. The thermocouple 
voltage is read by the computer and recorded as part of the 
data each time a measurement is made. 


Programs 

Four different variations of the Basic program were 
developed. (1) Magnitude, which measures the magnitude or 
the peak, at the sweep coil frequency, as a function of 
magnet current, from 0 to 1 A. (2) Cross Spectrum, which 

measures the real part of AxB, where A is the sweep coil 

drive signal and B is the transducer output signal. The 
field current is swept from 1 A to -1 A. (3) Quick, which 

measures the magnitude at five values of magnet current , 

from 0 to 0.5 A and calculates the average slope and (4) TC, 
which records the magnetization and temperature as a sample 
is warmed through the superconducting transition. The flow 
chart in figure 2 applies to all three programs. The 
programs are menu driven. The first section initializes the 
four instruments and sets the operating parameters. The 
operator selects the frequency and voltage for the sweep 
coils, the input voltage range for the analyzer, and the 
qain for the display Y-axis. In program 1 the option of 
field cooling is provided as well. Following the 
initialization the subroutine ''Acquire” is called which 
manages the actual acquisition of data. Once the data is 
taken, the subroutine "Store" is called and the data is 
stored on disk. The operator is prompted for a filename and 
may elect to not store the data at this point by doing a 
reset, in which case the data is lost. Once the data is 
stored on disk the menu provides options of plotting on the 
screen, printing on the screen, plotting on the plotter, or 
printing on the printer. Program listings are given in the 
Appendix . 
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Results 


Of the three types of piezoelectric transducers used, the 
PZT ceramic was found to give the most stable and 

reproducible results. This was largely due to the failure of 
the silver paint to properly bond the Kynar films to their 
copper supports. The response of the PZT reed to a white 
■ S ^ na l ■ a PPli e d to t ^ le swee P coils is shown in Figure 
3. This particular reed has a resonant frequency of 111 Hz. 
The peak at 10 Hz is due to vibrations of the isolation 
table on which the magnetometer was mounted. Figure 4 shows 
the response curve for an unsupported Kynar film. No 

resonance is seen, the 10 Hz noise peak is particularly 
large, followed by peaks at 60 and 180 Hz due to line noise. 
The response of this reed was nearly independent of 

frequency but was not nearly as good as that of the PZT reed 
near resonance. All of the following data were taken with 
the PZT reed. 


Initially, two types of curves were generated. The curves of 
magnetization vs magnetic field, and curves of 
susceptibility vs temperature, which were generated by 
taking average slopes of the magnetization curves, fitted to 
straight lines. Figures 5 and 6 show a set of magnetization 

fo 5 Er l Ba 2 Cu 3°x' the material which is refered to as 
Erl 23. The value of x is between 6.5 and 7 in the 
superconducting phase of this compound. In figure 5 the 
material is superconducting and diamagnetic while in figure 
6 it is normal and paramagnetic, due to the Er ion in the 
lattice. The average slope of each curve was taken to 
determine the change in magnetic susceptibility with 
temperature, the result is shown in figure 7. 


The degree of magnetic flux pinning can be determined by 
comparing magnetization curves for samples which have been 
f ?;? i “ < ? ooled ( Meissner effect ) and zero-field cooled ( 
shielding ) . If a sample is cooled below its superconducting 
critical temperature while in an applied magnetic field, 
the internal field is expelled, this is the Meissner effect. 
For type II materials, however, some magnetic flux can be 
pinned in the material and, for small fields, will not be 
expelled. This results in a positive contribution to the 
magnetization. Thus a difference in slope between field- 
cooled and zero-field cooled data on the same sample will 
give an indication of the extent of the flux pinning. 

Figure 8 shows the magnetization vs magnetic field curves 
fS- 97 a ® a ^ n P le of Y123. The curves represent field-cooled 
0e ) . and zero-field cooled (H=0) conditions. The 
difference m slope between the two curves indicates the 



Frequency Response of Reed # 3 to Noise 7/11/89 AVERAGE 




Figure 3. Frequency Response of the PZT Reed. 
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Figure 7. Susceptibility Curve for Erl23. 
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extent of the flux pinning. Figure 9 shows the same 
experiment for Y123 which has been doped with Ag. Here the 
slopes of the two curves are practically identical, 
indicating little flux pinning in the doped material. 


The program "Cross Spectrum" has the capability of 
distinguishing between positive and negative magnetization 
and can be used to generate hysteresis curves for magnetic 
materials. Figure 10 shows a hysteresis curve for a sample 
of Eul23+Ag which was zero-field cooled. The increase m 
(negative) magnetization on going once around the loop xs a 
measure of the amount of trapped flux. The same sample is 
shown in figure 11 after field-cooling in a 66.4 Oe magnetic 
field. The curve is now closed, indicating that flux was 
trapped upon cooling . 
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Figure 9. Magnetization Curve for Field-Cooled Y123+Ag 


XII-16 










Field current CR) 



Conclusions and Recommendations 

The alternating gradient magnetometer can be a valuable tool 
for thf characterization of the magnetic f 

materials. It should be P^ rt ^ a {%i ly Bgn sltivlty and 
simerconductivity research because of itB sensitivity »■ 

of use With the current data acquisition system, data 
may be £ taken "ore quickly than before and the data analysis 
T! Y areatly simplified. The data which is acquired will 
become even more useful if the instrument is calibrated to 
measure magnetization in absolute units. This should not be 
St«Icult to do using a standard material, such as one of 
the rare-earth salts* 

The chief difficulty in using this instrument 18 
precise temperature 7 control . The temperature of the sample 
changes with the liquid nitrogen level in the dewar^ As the 
liauid evaporates the sample slowly warms. It wouia oe 
convenient to be able to take a series of magnetization 

measurements at a fixed temperature This is not possible 

with the current instrument although the rate of temper 
chanai can be kept quite low, so that the temperature does 
nSr change more than about 0.1 K during the course of a 
measurement . 

The computer program can be improved in several places . The 
proaram called Magnitude could be rewritten so that it would 
take a series of measurements, say one ^magnetization Q ^ilk 
each degree of temperature change, and save them to dis 

automatically. The program called Cross A^toTt aTtell 
have the temperature measuring routine added to it as w • 
Finally, some way needs to be devised to keep track of any 
variation in the resonant frequency of the reed and adjust 
the function generator accordingly. 
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1 0 
20 
21 


I X X- -X X X X X * X -X X X X- X X X X X X- X X X X X- X X * X- X X X X X X- X X X X X X X X X X X X X X X X X' X X X X X 

PROGRAM : MAGNITUDE 


30 

! PLOTS THE 

MAGNITUDE OF THE SIGNAL ON CHANNEL B 

! 

40 

! FROM 

ZERO TO 1/2 AMP MAGNET CURRENT 

1 

50 

* USING 

THE CF-920 FFT ANALYSIS SYSTEM 

i 

60 

j 

JUNE 8, 1989 

! 

70 

i 


i 

B0 

! MODIFIED 

JULY 7, 19B9 TO USE THE PZT TRANSDUCER 

i 

81 

! MODIFIED 

JULY 17, 1989 TO MONITOR TEMPERATURE 

) 

B2 

! 


1 

90 

! xxxxxxwxxxxx*******#***##*****#* **#**#*#***#### ********* ! 

1 00 

i 



1 1 0 

! Nunber of Averages is in Line 580 


120 

j 



130 

OPTION BASE 1 



140 

CONTROL t , 1 2 j 1 



150 

DIM Magnitude(250).Fieldcurrent ( 250 ) , Tenpera t ore( 250 ) 


160 

PLOTTER IS 3," 

INTERNAL" 


170 

GCLEAR 



180 

! 



190 

! ****** ******** 


200 

! 



21 0 

PRINT CHR $(12) 



220 

PRINT ” (1) 

Acquire nagne t i za t i on data" 


230 

PRINT 



240 

PRINT *’ (2) 

Load Magnetization data froM disc" 


250 

PRINT 



260 

PRINT M < 3 ) 

Plot Magnetization data on screen” 


270 

PRINT 



280 

PRINT '* (4) 

Print Magnetization data on screen” 


290 

PRINT 



300 

PRINT " (5) 

Plot Magnetization data on plotter” 


310 

PRINT 



32 0 

PRINT " (6) 

Print Magnetization data on printer” 


33 0 

INPUT Selection* 


34 0 

IF Selections 

"1” THEN GOSUB Initialize 


350 

IF Selection$= 

”2” THEN GOSUB Load 


360 

IF Selections 

"Z" THEN GOSUB Screenplot 


370 

IF Selections^ 

”4” THEN GOSUB Screenprint 


380 

IF SelectionS= 

”5" THEN GOSUB Plotterplot 


390 

IF Selection$= 

”6” THEN GOSUB Printerprint 


400 

GOTO 210 

! Clear Screen and Show 

Men 


^110 Initialize: ! Setup Display Window 

420 PLOTTER IS 3 , " INTERNAL " 

GINIT 
GCLEAR 

VIEWPORT 0,100,0,100 
WINDOW 0,1 ,0,3 
GRID ,25,1 
GRAPHICS ON 
FRAME 


430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 


I 

i xxxxxxxxxxxxxx******* SET UP FOR CF920 **#**#*#«##«#*******#***** 


1 


INPUT "Channel B Ranqe <l-15> H ,Br$ 

IF RrS=”” THEN Br$*"10” 

OUTPUT 710;"AMS3" ! Spectrun peak averaging 
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580 

OUTPUT 71 0 ; "ANS5" 

! Take 5 squared averages 


590 

OUTPUT 71 0 ; " TR E H 

! Repeated trigger 


600 

OUTPUT 7 1 0 j " T I A " 

! Channel a trigger 


6 1 0 

OUTPUT 71 0 ; " AAS2 " 

! Channel A input voltage range 


620 

OUTPUT 710; “BAS ”&BrS 

! Channel b input voltage range 


630 

OUTPUT 710;"FRS8 M 

! Channel A&B input frequency range 


640 

OUTPUT 710 ; "BSP " 

! Display B spectrum 


650 

OUTPUT 710) "YLI" 

! Selects linear scaling 


660 

OUTPUT 71 0; “SON- 

! Search function on 


670 

OUTPUT 71 0)"SPS401 ’* 

1 Sets search marker for peak point 


660 



! 

690 

! ttttttff* ****#### ***** CF9 

20 SETUP COMPLETE ************************* 

700 



\ 

710 



! 

72 0 



| 

730 

! * * * x x * x x- * * x x * * * # * *SET U P 

FOR THE 228 CURRENT SOURCE**************** *** 

740 



! 

75 0 

CLEAR 711 



760 

OUTPUT 711;”V25" 

! Voltage limits to 

25 volts 

770 

INPUT "Initial Current <0 

-1> M , Current 1 Used for Field 

Cooling 

780 

IF Currents THEN GOTO 810 


790 

OUTPUT 71 1 ; " I "4VAL* (Current ) 


800 

GOTO 820 



BIO 

OUTPUT 71 1 ) "10 " 

1 Initializes current to 0 

amperes 

820 

OUTPUT 71 1 ; "FIX" 

! Execute 


830 

INPUT "Press a Key to Continue". K$ 


840 

1 x -x xxx x- x x x x x x SETUP FOR 

THE 229 CURRENT SOURCE COMPLETE ****x***** 

850 



! 

B60 

!************ SETUP FOR 

HP3457A TO READ TEMPERATURE ***************** 

870 

! USING A COPPER-CONSTANTAN THERMOCOUPLE 

1 

eso 

! 



B90 i 

INTEGER I 



900 

ASSIGN gTenp TO 724 



910 

OUTPUT 9Temp j "TRIG SYN M 



920 

OUTPUT 8Temp; "NDIG 4" 

! Four Digit Accuracy 


930 

OUTPUT BTenp j "FIXEDZ OFF" 



940 

OUTPUT @ T e mp ; "DCV" 

i Read DC Volts 


950 

OUTPUT PTenp ; "AZERO 1" 

! Automatic Zero 


960 

OUTPUT PTemp j "NPLC 1 0 " 

! Ten Power Line Cycle Delay 


970 

j 



980 

! x-K«***tt*tt****tt***ft*ft* TEMPERATURE SETUP COMPLETE ***#**#****#*#****** 

990 

! 



1000 

! x xxxx -x -x- ***** xxx x x-x-x SETUP FOR THE 270 * * x x x x x x * xxxx x x x ******** x x 

1010 



i 

1 020 

Rf $="1 00 ” 

1 Resonant frequency of PZT 


1 030 

Swv$=" 4 , 40" 

! Sweep voltaqe of little coils 


1 040 

PRINT CHR $ (12) 

! Clears screen 


1 050 

CLEAR 709 



1 060 

INPUT "Resonant frequency 

",Trf* 


1 070 

IF Tr f $= " " THEN 



1 080 

Rf$=Rf$ 

! Default value is 74.4 Hz 


1090 

ELSE 



1100 

Rf $=Trf $ 



1110 

END IF 



1 120 

INPUT "Sweep v o 1 1 age M , Ysu vS 


1130 

IF T swv$= " " THEN 



1 140 

Suv$=Swv$ 

\ Default value is 4.4 volts 


1150 

ELSE 



1160 

Su v$=T swv$ 



1170 

END IF 
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1180 
1190 
120 0 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1 390 
1400 
1 41 0 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
151 0 
1520 
1530 
1540 
155 0 


Freq=VAL(Rf$) 
Rf$="F"6Rf$ 
OUTPUT 70 9 ; Rf $ 
OUTPUT 709; "I" 
Suv$=" A ”6Swv* 
OUTPUT 709; Swvf 
OUTPUT 709; "I" 
OUTPUT 709; "PI- 
OUTPUT 709; "I" 
OUTPUT 709; "TSI" 
OUTPUT 7 09 ; "F " 

! 


Adds a F for frequency to variable 
Sets frequency of Uavetek 
Execute 

Adds an A for amplitude to variable 

Sets voltage for Uavetek 

Exec u te 

Output on 

Execute 

Displays frequency on front of 270 


!************* ENABLE CF920 DELTA SEARCH FUNCTION 


OUTPUT 710; "SON " 

Hf = (Freq ) *2+5 ! 

Lf = < Freq ) *2-5 I 

Lf i=VAL$ ( Lf ) 

Hf $=VAL$ ( Hf > 

OUTPUT 71 0;-SPS-6Lf$ 
OUTPUT 71 0; "DON- 
OUTPUT 710; "DCS - 
OUTPUT 71 0 ; "SPS'^Hf* 
INPUT "Y-AXIS GAIN <l-5> 
FOR G = 1 TO Yag 
OUTPUT 710; " YGU M 
NEXT G 

OUTPUT 710 ; "PAQ M 


I SEARCH ON 

UPPER LIMIT OF SEARCH 

LOWER LIMIT OF SEARCH 


1 SET LOWER LIMIT 
I DELTA FUNCTION ON 
! SET THE DELTA CURSOR 
! SET UPPER LIMIT 
,Yag 


! PARTIAL OVERALL ON 


****** *»***«•****• j 


! 

! * * * * * * * x- x- x -x- * * SETUP FOR THE 270 UAVETEK COMPLETE **************** 

\ 

\ 

j 

GOSUB Acquire ! Take data 

! ***«*«***»***********««ft*ff»**ft***«»*tf**********»***tt«**»«»*«*«««*»*ftft 

GOSUB Store i Write the data to disc 

RETURN 1 Back to the Menu 

STOP 


1560 
1570 
1580 
1590 
1600 
1 61 0 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
171 0 
1720 
1730 
1740 
1750 
1760 
1770 


Acquire: ! Take data ************************************************* 

1 = 0 


FOR Fieldi=0 TO ,5 STEP ,1 ! 0 to 0.5 Amps in 5 steps 

IF FieldiU.E-9 AND F ield i >-l . E-9 THEN Fieldi=0 
1 = 1 + 1 


F i e Id i " I " 6V ALS ( F ie 1 d i ) ! Adds an I to the field current vairable 


OUTPUT 71 1 ; Fieldi$ 
OUTPUT 711; "FIX- 
WAIT 3 

OUTPUT 710; H LYS- 
ENTER 71 0 ; Vo 1 ts 
Magnitude(I)=»Volts 


! Sends value of current to current source 
! Execute 

! Wait for system to settle 
! Selects Display Y search value 
Reads display Y search value ; volts 


Fieldcurrent<I)«Fieldi 

GOSUB Plot \ Plot data on the screen 

NEXT Fie Id i 

FOR Fieldi= . 5 TG 0 STEP -.1 I .5 to 0 amps 

IF Fieldi<l.E-9 AND Field i >-l . E-9 THEN Fieldi=0 
1 = 1 + 1 


Field i$=" I "6VAL$< Field i) ! 
OUTPUT 711 ; Fieldit ! 
OUTPUT 71 1 ; "Fl X - ! 
WAIT 3 ! 


Adds an I to the field current vairable 
Sends value of current to current source 
Execute 

Wait for system to settle 
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1780 
1790 
180 0 
1810 
1820 
1830 
1840 
1650 
1860 
1870 
t B8 0 
1890 
1900 
191 0 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
200 0 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
21 1 0 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
221 0 
2220 
2230 
2240 
2250 
226 0 
2270 
2280 
2290 
230 0 
2310 
2320 
2330 
2340 
2350 
2360 
2370 


OUTPUT 710) "LYS" ! Selects display Y search value 

ENTER 71Q;Volts ! Reads display Y search value 
Magnitude* I )=Volts 
Fieldcurrent ( I)=Fieldi 


GOSUB Plot 
NEXT Fieldi 
ENTER BTemp ;Voltage 
Volts=Voltage*1000 
GOSUB Temperature 
OUTPUT 71 1 ) "FOX" 
CLEAR 709 
OUTPUT 710) " P AF * 
OUTPUT 710; “DOF N 
OUTPUT 71 0 ; H SOF" 
CLEAR 710 
RETURN 


! Show data point on screen 

! Read the Temperature 
! Convert to Millivolts 
! Calculate the Temperature 
! Put Current Source on Standby 

! DISABLE PARTIAL OVERALL 
! DISABLE DELTA SEARCH 
! DISABLE SEARCH FUNCTION 

\ From ACQUIRE 


Plot: ! Plot data points on the screen **************************** 

LORG 5 ! Label origin at center 

MOVE FieldcurrenT(I)>Magnitude(I) 

LABEL ! Plot using asterisks 

PRINT I fieldcurrent ( I ) ,Magni tude< I ) 

RETURN 


Store: ! Store the data in a disc file ****************************** 


N=I 


REDIM Magnitude(N) fieldcurrent(N) 

GRAPHICS OFF 
ALPHA ON 

PRINT CHR $ C 12) » Clear screen 

INPUT “File name ?"file* 

CREATE BDAT File*, INT ( N*8/256 ) +2 
ASSIGN ePath TO File* 

OUTPUT PPath)N,Temp,Magnitude(#> fieldcurrent**) 
ASSIGN 8Path TO * 

RETURN ! From STORE 


Load: ! 

PRINT CHR $ ( 1 2 ) 

Control ler *1 
CAT ) NAMES 

INPUT "Enter file name to 
IF File**"" THEN 2240 
ASSIGN BTrack TO File* 

ENTER PTra c k ) N 
ENTER PTrack)Temp 

REDIM Magn i t ude<N ) fieldcurrent (N) 

ENTER BTrack; Magnitude**) fieldcurrent (*> 
ASSIGN BTrack TO * ! Close the file 

RETURN ! From LOAD 


! Clear screen 

! Display catalog filenames 
be loaded from above list, "file* 
! Return if no entry 

! Number of data points 


Sc r eenpr in t s ! Print the data on the screen 
IF Con tr o 1 ler = 0 THEN 


DISP M LOAD A FILE FIRST!! I" 

GOTO 2480 
END IF 

PRINT "Filename « "file* 

PRINT “Temper a t ur e » M ;Temp 
FOR 1=1 TO N 

PRINT I,Magnitude(I)fieldcurrent<I) 

NEXT I 

FOR K = 1 TO 4 

Mag(K)=Magnitude(KM )-Magnitude(l ) ! Calculate slope 

NEXT K 
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2380 

Su «=0 


2370 

FOR K = 1 TO 4 


2400 

Sun=Sun+Mag(K ) 


2410 

NEXT K 


2420 

SI ope=Su*/ . 5 

! Average slope of first 5 poii 

2430 

PRINT 


2440 

PRINT "Initial Slope = 

" j S 1 o p e 

2450 

DISP "Hit CONTINUE to 

continue progran." 

2460 

PAUSE 


2470 

Control ler=0 


2480 

RETURN 

! -Fron SCREENPRINT 

2490 

1 


2500 

Pr in terpr int : ! Send 

data to the printer 

2510 

PRINTER IS 701 


2520 

GOSUB Screenprint 


2530 

PRINTER IS 1 


2540 

RETURN 

! Fron PRINTERPRINT 

2550 

1 


2560 

Scr eenp 1 o t : 1 Plot 

the full graph on the screen 

2570 

PRINT CHR $ ( 1 2) 

! Clear screen 

2580 

GCLEAR 

1 Clear graphic screen 

2590 

GRAPHICS OFF 


2600 

ALPHA ON 


2610 

GOSUB Sort 

i Find nin and wax values 

2620 

GRAPHICS ON 


2630 

ALPHA OFF 


2640 

Yinta=Yint-<Maxy-hiny >/75 

2650 

Xinta=Xint-(Maxx-Minx >/75 

2660 

VIEWPORT 20,100,12,92 


2670 

FRAME 


2680 

WINDOW hinx,Haxx ,Miny , 

Maxy 

2690 

AXES Xtic , Ytic ,Minx ,hiny , 1 , 1 ,2 

2700 

PEN 1 


271 0 

LOR G 6 

! Label origin above character 

2720 

GRID Xtic, Ytic 

! Draw a grid 

2730 

CLIP OFF 


2740 

CSIZE 3, .5 

* Character size/aspect ratio 

2750 

FOR I = Min x TO ha x x+X t ic /t 0 0 STEP Xtic 

2760 

MOVE I , Min y 

t X-axis nunbers 

2770 

IF ABS( I ) < 1 . E-15 THEN 

2780 

J=Q 


2790 

LABEL J 


28 U 0 

GOTO 2830 


281 0 

END IF 


2820 

LABEL I 


2830 

NEXT I 


2840 

LORG 7 

! Label origin lower right 

2850 

FOR I=Miny TO Ma x y +X t ic/1 0 0 STEP Ytic 

2960 

MOVE Minx, I 

t Y-axis nimbers 

2870 

IF ABS(1)<1 .E-15 ' 

THEN 

2880 

J=0 


2890 

LABEL J 


2900 

GOTO 2930 


2910 

END IF 


2920 

LABEL I 


2930 

NEXT I 


2940 

LORG 5 

! Label origin to center 

2950 

PEN 2 

! Pen for data points 

2960 

FOR 1=1 TO N 


2970 

MOVE Fieldcurrent(I) , Magnitude(I) 
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2980 
2990 
30 0 0 
3010 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
, X t i c 
3380 
3390 
3400 
3410 


LABEL "o" ! Character for points 

NEXT I 
PEN 1 

VIEWPORT 20,100,15,95 
WINDOW 20,100,15,95 
CLIP OFF 

MOVE 60,5 * Position of X-axis label 

LORG 5 

LDIR 0 1 Horizontal 

LABEL "Field current <A)“ ! X-axis label 

MOVE 5,55 ! Position of Y-axis label 

DEG 1 Use degrees 

LDIR 90 ! Orient vertically 

LORG 5 ! Center 


LABEL "Magnetization (mV) 
LDIR 0 
MOVE 115,75 
LABEL Title* 

LABEL "Tenp, - "jTemp 
LABEL "Filename * ";File* 
LABEL TIME* ( TIME DATE ) 
LABEL DATE* ( T IMED ATE ) 


I Y-axis label 
t Hor i z on t a 1 
I Position for sidebar 
! Sidebar header 


I Time 
! Date 


RETURN ! from SCREENPLOT 

I 

Sort : ! Get max and nin 

Minx=MIN(Fieldcurrent(*) ) 

Maxx=MAX<Fieldcurrent(*)) 

Mmy = hIN(Magnitude(*) ) 

Maxy=MAX(Magnitude(*) ) 

PRINT "F ieldc urren t ranges f r o m " \ Mi n x ; " t o " ; Ma x x 
PRINT 

PRINT "Magnetization ranges f r om " j Mi n y \ " t o " j Ma x y 
PRINT 

INPUT "Input ( Min x , Ma x x , Min y , Ma x y ) to be p 1 o t t ed " , Min x , Ma x x , Mi n y , Ma x y 
PRINT 

PRINT "Fiel dc urren t plot ranges f rom" ; Min x ; " t o " ;Max x 
PRINT 

PRINT "Magnetization ranges plot f r o m " j Mi n y j " t o " ; Ma x y 
PRINT 

INPUT "Input tick mark spacing and axes intercepts (Xt ic , Ytic ,Xint , Yin t ) " 
, Y Tic , Xin t , Yin t 
PRINT 

INPUT "Title for graph ", T i 1 1 e* ! Sidebar title 

RETURN I From SORT 

i 


3420 P 1 o t t erp 1 o t : \ Send plot data to plotter 

3430 PLOTTER IS 705,"HPGL" 

3440 GOSUB Screenplot 

3450 RETURN ! From PLOTTERPLOT 

3460 Temperature: I Calculate the Temperature from Thermocouple Voltage 

3470 LET V=Vo 1 ts 

3480 IF V > =-5 ,60 2 AND V<-5,43? THEN 
3490 T emp = 73+61 . 1 09* ( V+5 . 603) 

3500 END IF 

3510 IF V>=-5 , 439 AND VC-5.261 THEN 
3520 Temp=83+58 . 602# (V+5 . 439) 

3530 END IF 

3540 IF V >— 5 . 261 AND V<-5.069 THEN 
355 0 Temp = 93 + 54. 1 59* < V+5. 261 ) 

3560 END IF 
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3570 
3580 
359 0 
3600 
3610 
3620 
3630 
3640 
3650 
3652 
3660 
3670 


IF V >-5 . 069 AND U<-4.865 THEN 
TenpM 03+50 . 596 * < V+5 . 069 > 

END IF 

jF V>-4.865 AND V<-4.468 THEN 

Tenp=113+47.507*(V+4. 865 ) 

END IF 

IF y>*- 4.468 THEN 

Tenp=*V 

END IF 

PRINT Temp 

RETURN 

END 


I Just report the voltage 


! From Tenperature 
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